Fermi Surface Reconstruction by Dynamic Magnetic Fluctuations 
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We demonstrate that nearly critical quantum magnetic fluctuations in strongly correlated electron 
systems can change the Fermi surface topology and also lead to spin charge separation (SCS) in 
two dimensions. To demonstrate these effects we consider a small number of holes injected into 
the bilayer antiferromagnet. The system has a quantum critical point (QCP) which separates 
magnetically ordered and disordered phases. We demonstrate that in the physically interesting 
regime there is a magnetically driven Lifshitz point (LP) inside the magnetically disordered phase. 
At the LP the topology of the hole Fermi surface is changed. We also demonstrate that in this 
regime the hole spin and charge necessarily separate when approaching the QCP. The considered 
model sheds light on generic problems concerning the physics of the cuprates. 

PACS numbers: 74.40. Kb, 74.72.Gh, 75.10.Jm, 75.50.Ee 
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It is well known that static magnetic order in a conduc- 
tor can influence the Fermi surface (FS) topology due to 
the electron/hole diffraction from the ordered moments 
(Fermi surface reconstruction) . It is also well known that 
spin and charge are separate in one dimensional systems 
[3, 0] . In the present work we address two generic prob- 
lems - (1) Can dynamic magnetic fluctuations drive a 
change in the FS topology; (2) Is it possible to separate 
spin and charge in a two-dimensional (2D) system, and 
if so what is the meaning of the separation. We give 
positive answers to both questions and demonstrate that 
these two problems are remarkably related. To address 
the generic problems we consider a specific model of a 
small number of holes injected into the bilayer antifer- 
romagnet (AF) with magnetic fluctuations driven by the 
interlayer coupling. We show that indeed purely dynamic 
short range AF correlations in the absence of a static AF 
order can cause a LP in which the topology of the FS 
is changed. A similar model has been analyzed previ- 
ously by Vojta and Becker Q where they also observed 
a LP. However in their case the LP was always in the 
AF ordered phase and therefore the central issue of the 
magnetic fluctuation driven LP has not been addressed. 
We also demonstrate that when the LP is in the disor- 
dered phase the hole spin and charge necessarily separate 
when approaching the magnetic QCP. The possibility of 
spin-charge separation (SCS) in 2D has been discussed 
previously within the context of the slave-boson method. 
The method applied to the t — J model implies the sepa- 
ration ad hoc [J] . In the present analysis we do not make 
any ad hoc assumptions. The precise meaning of the sep- 
aration following from our analysis is different from that 
of the slave-boson method. 

Our interest in the problems is motivated by the 
cuprates, which, in our opinion, manifest both the LP 
and SCS. Below we explain the connection with cuprates. 
A reader who is not interested in cuprates can go directly 
to Eq.(l) where we start the analysis. 



Lying at the center of the debate of high-T c supercon- 
ductivity is whether it originates from a Fermi liquid or 
from a Mott insulator. The angle-resolved photocmis- 
sion spectroscopy (ARPES) indicates that the transition 
from a small to a large FS occurs in the hole doping range 
0.1 < x < 0.15 Existence of hole pockets is con- 

sistent with the picture of dilute holes dressed by spin 
fluctuations, based on doping a Mott insulator [8[. On 
the other hand, ARPES in the optimally to overdoped 
cuprates show a large FS as expected by the Fermi liquid 
approach. This implies that there is a topological Lifshitz 
point [9] (LP) in the doping range 0.1 < x < 0.15, where 
the Fermi surface changes from small to large. A phe- 
nomenological description of the LP based on the Fermi 
liquid picture was suggested in Ref. [l(|. DMFT calcula- 



tions 



llj with the Hubbard model also indicate the LP. 



Magnetic quantum oscillation (MQO) in underdoped 
YBa2Cu30 y [12( have revealed a small FS pocket, in con- 
trast to the large FS observed on the overdoped side [l]| ■ 
This again indicates the existence of a topological LP. 
The MQO measurements were performed in very strong 
magnetic fields, up to 80T. Therefore, a possible point of 
view is that the field induces a static magnetic structure 
and the structure causes the small Fermi surface recon- 
struction [l4[. However, the small FS was observed in 
MQO up to 12% doping [THJ] , and it is unlikely that even 
an 80T field can generate a static AF order at such high 
doping. On the other hand the short range dynamic AF 
correlations always exist in the cuprates. Moreover, re- 
cent RIXS measurements (l6j demonstrate remarkably 
that such correlations are practically doping indepen- 
dent, from Mott insulator to optimal doping. Based on 
this data one can conjecture that the cuprates are always 
close to magnetic criticality. This motivates us to study if 
the LP can be driven by short range, purely dynamic AF 
correlations. We consider a bilayer model for the sake of 
performing a controlled calculation. However, we believe 
that conceptually our conclusions are equally applicable 
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to both single and multi-layer cupratcs. 

Now we turn to the discussion of SCS and its relation to 
magnetic quantum criticality. Optimally and overdoped 
cuprates do not have any static magnetic order. On the 
other hand, the underdoped cuprates possess a static in- 
commensurate magnetic order at zero temperature. A 
magnetic QCP separating these two regions was pre- 
dicted in [17|. In La2_ :E Sr a ;Cu04 the QCP is smeared out 
because of disorder. However, in YBa 2 Cu 3 O y the QCP 
is located experimentally at doping x ~ 0.09 (y « 6.47) 
[l8l - l20j . In the magnetically ordered phase the hole does 
not carry usual spin, instead it carries only a psudospin 
which marks the sublattice. The pseudospin interacts 
with a magnetic field in a very unusual way [HI [22| and 
this is the mea ning of the partial SCS in the magnetically 
ordered phase [17|. On the other side of the QCP there 
is no static magnetic order, and hence spin and charge 
are united. In the present work we analyse the process 
of SCS at the QCP. The model considered here has only 
commensurate magnetic ordering, so we put aside incom- 
mensurability in the cuprates. 

We consider the t — t' — t" — J model defined by the 
following Hamiltonian on the bilayer square lattice 

~ ^UAcLcj" + c )a c *<y) > (!) 



where c i<T is the creation operator of an electron with spin 
a =|, I at site i on the top plane, = <t ^ u c iv , and 
tij = {t,t',t"} is the hopping integral between nearest, 
next-nearest, and next- next nearest neighbour sites re- 
spectively. The superscripts (1), (2) indicate the layers. 
Hereafter we set J = 1. A no-double-occupancy con- 
straint is imposed. It is well known that without holes 
(half- filling) the model has an 0(3) magnetic QCP at 
J_L ~ 2.5 [23J separating the magnetically disordered and 
the AF ordered phases. The hopping integrals t, t', t" re- 
sult in charge dynamics if a hole is injected into the sys- 
tem. The longer range hopping integrals t' , t" are crucial 
as we will explain later. Note that we consider the zero 
temperature case, so the magnetic ordering in the AF 
phase is consistent with the Mcrmin- Wagner theorem. 

We focus on small doping, x <C 1, such that it does 
not influence the magnetic fluctuations. Magnetic fluctu- 
ations and the QCP are driven by the interlayer coupling 
J_l • The holes fill the rigid band formed by the magnetic 
quantum fluctuations. To address the problems formu- 
lated above it is sufficient to calculate the single hole 
Green's function. Certainly at sufficiently high concen- 
tration of holes they start to influence the magnetic fluc- 
tuations and hence the rigid band approach fails. How- 
ever, we do not need to go to such high concentrations 
to draw our conclusions. Such an approach is only pos- 
sible because the magnetic dynamics are driven by Jj_ 



and arc independent of the hole concentration. This is a 
significant simplification compared to the t — .//Hubbard 
model, where doping is the only "handle" . 

We find that close to the 0(3) QCP, the hole dis- 
persion has minima at k = (±-|,±-|). This results in 
small hole pockets similar to that in the cuprates at small 
doping. The pockets are formed due to strong in-plane 
AF correlations which diminish the nearest site hopping 
t. Upon increasing J± the in-plane AF correlations are 
reduced and the dispersion minima gradually shift and 
reach k = (±7r, ±7r) at some value J± p ■ This is the 
position of topological LP at x — ► 0. The in-plane AF 
correlations are diminished at J±_ — > oo and hence the 
hole dispersion is almost like that in the normal Fermi 
liquid, but a factor of two reduced due to J±. 
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2t lk + 2t' lk + 2t" lk 



(2) 



■ [cos k x + cos k y ] 



7 fe = cosk x cosk y , 7 fe = — [cos(2fca;) + cos(2fcj,)] 

At J± < J^ p there are four FS pockets and at J± > J\ p 
there is only one pocket centered at (7r,7r) according to 
([2]). It is important to make it clear that as soon as we 
consider only one hole, we observe just a crossover from 
one shape of the dispersion to another. However, at a 
finite concentration of holes this results in a sharp LP 
with MQO frequency jump by a factor of four. 

Magnetic excitations in the magnetically disordered 
phase are triplons. To describe the triplons we employ 
the spin-bond operator mean field technique This 
technique gives a QCP at Jj_ « 2.31 which is reasonably 
close to the value 2.525 known from Quantum Monte 
Carlo [23| . All necessary Eqs. describing the triplon 



dynamics have been derived in Rcfs. [24|. One can cer- 
tainly employ a more accurate Brueckner technique [25| . 
However, this technique is more involved while the bond 
operator mean field approach has sufficient accuracy for 
our purposes and we chose it for simplicity. 

To describe the hole dressed by triplons we use the self- 
consistent Born approximation (SCBA) which disregards 
vertex corrections. This approximation has been widely 
used to study hole dynamics in the AF background [8|. 
In the AF background the single loop vertex correction 
is zero due to kinematic constraints and hence the SCBA 
is very accurate. In the present case of the magnetically 
disordered background the single loop vertex correction 
is nonzero. However, the correction is suppressed by the 
parameter 1 /N, where TV = 3 is the number of the triplon 
components [26j. To confirm the accuracy of the SCBA 
we compare results with that of numerically exact dimcr 
series expansions. Note that here we are working in terms 
of the true spin of the hole, while in the case of the AF 
background one has to work in terms of pseudospin. 

The SCBA results in the following Dyson's equation 
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FIG. 1: (color online) The hole spectral functions at t — 
0.5, t' = t" = 0, for J± = 2.31 (a) and J ± = 3 (6). Here 
we show k = (-|, ^) (black, solid), (7r,7r) (red, dashed), and 
(0,0) (blue, dotted). There is a significant distant incoherent 
part, the quasiparticle residues for J± — 3 are Z^K t s.) = 0.66, 
^(tt.tt) = 0.88, Z(o,o) = 0.32. 



for hole Green's function 

G(k, e)=(e- - 3 £ <£_ q , q G(k - q, e - fi fl ) + iO 



where the factor 3 comes from the three different po- 
larizations of the intermediate triplon. The bare hole 



dispersion ef 1 ' is given by 1^. 

1S 9k,q = -77}7 [UqT 



T p = [2i 7p 



The hole-triplon vertex 
^ where 



V q T k + q ] - ^lq{u 

and v. 



and 



are the triplon 



Bogoliubov coefficients. These coefficients as well as the 
triplon energy Q q were calculated in [24[. Note, that on 
approaching the QCP the vertex diverges at q — > (tt, tt) 
and this leads to the SCS as we discuss below. A similar 
Dyson's equation was used in Ref. 
ferent vertex gk,q- We believe that the vertex in 



271, but with a dif- 



27j is 



wrong. We solve Eq. §3§ numerically on a 128 x 128 clus- 
ter with energy resolution Ae = 0.02. The quasiparticle 
dispersion is given by the position of the (^-function peak 
in the hole spectral function A(k, e) = —-Im[G(k, e)]. 

We start from the case of small hopping, t = 0.5, 
t' = t" = 0. Wc plot in FigJT] the spectral functions 
at values J± = 2.31 and J± = 3. Note that J± = 2.31 is 
exactly the position of the QCP obtained from the mean 
field triplon analysis. Spectra in FigOJa) do not show 
any quasiparticle peaks, instead there are only power 
cuts. This is similar to the Green's function of an immo- 
bile magnetic impurity at the QCP [26, 28 1. Remarkably 
the hole mobility does not influence this behaviour. The 
power cuts imply that the spin is distributed around the 
hole in a diverging cloud indicating SCS at the QCP, see 
discussion below. On the other hand spectra in Fig. [TJb) 
show quasiparticle peaks separated by the triplon gap A 
from the incoherent spectra. Figure [ljb) shows the dis- 
persion minimum at k = (n, it), while the cut position in 
Fig. [TJa) is practically the same for all momenta. Hence, 
we conclude that the position of the LP in this case co- 
incides with that of the QCP. 

We also perform dimer series expansion calculations 
for t = 0.5, t' = t" = 0, and compare the results with 
the SCBA. The series expansion method allows one to 
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FIG. 2: (color online) The hole spectral function at t — 3.1, 
t' = t" = 0, for J± = 2.31 (a) and J± = 3 (6). Values of the 
momentum are k = (£, -|) (black, solid), (ir,ir) (red, dashed), 
and (7r,0) (blue, dotted). The quasiparticle residues for J± = 
3 are Z^2l : 2l) — 0.31, Z^^) — 0.80, Z(tt,o) = 0.35. The inset in 
Fig. (a) shows spectral functions in the broader energy range 
and the inset in Fig.(b) shows the map of the hole dispersion. 
The dark region is minimum of the dispersion. 



determine only the quasiparticle dispersion. Naturally 
the method does not converge close to the QCP, as there 
are no quasiparticles there. However, at Jj_ = 3 the 
method works well and agreement between the SCBA 
and series is good, for example the SCBA band width is 



=(0,0) 



0.40 and the scries band width is 0.41. 



In the strong coupling limit, t = 3.1, we rely on 
the SCBA since the series expansion does not converge. 
Spectral functions for t = 3.1, t' = t" = arc shown 
in Fig. [2] for J± = 2.31 and J± = 3.00. In this case 
there is no LP in the disordered phase, as the bottom 
of the band is always at k = (tt,tt). Hence the LP is 
inside the magnetically ordered phase in agreement with 
At the band bottom there are well defined quasipar- 
ticles even at the QCP. This indicates that there is no 
SCS. The spectra at k = (f , |) and fc = (tt, 0) have cuts 
at the QCP. However, these are the high energy states 
which are irrelevant at small doping. 

The last and the most important set of parameters, 
t = 3.1, t' = —0.8, t" = 0.7, roughly corresponds to the 
parameters of the cupratcs. Although wc do not intend 
to simulate the cuprates quantitatively, we choose the 
parameters relevant to cuprates such that the interplay 
of various mechanisms can be compared on a similar en- 
ergy scale. The spectral functions are shown in Fig. [3] 
for several values of J±. The dispersion maps shown in 
the insets clearly demonstrate that the LP is located at 
Ji m 3 within the magnetically disordered phase. This 
topological transition is caused by fully dynamic antifcr- 
romagnctic correlations. The demonstration of the pos- 
sibility of the fully dynamic scenario is the first major 
conclusion of the present work. Our analysis assumes 
small doping, practically we need x < 0.1 when the FS 
built on maps Fig. 3(b) and Fig. 3(c) are topologically 
different. Why the longer range hoppings have qualita- 
tively changed the situation? At t' = t" = the LP is in 
the ordered phase but already rather close to the QCP. 
A hopping t" > pushes the bare dispersion ([2|) at the 
nodal point k — (ir/2, 7r/2) down helping magnetic fluc- 
tuations to form a small pocket. A pretty small positive 
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FIG. 3: (color online) The hole spectral function at t — 3.1, 
t' = -0.8, t" = 0.7 for J x = 2.31(a), 2.50(b), 3.00(c) and 
4.00(d). Values of momentum are k = (§, ■§) (black, solid), 
(tt,tv) (red, dashed), and (7r,0) (blue, dotted). The quasipar- 
ticle residues for J± — 3 are Z^l^ = 0.29, Z^^) = 0.77, 
Z(ir,o) = 0.20; The inset in Fig. (a) shows spectral functions in 
the broader energy range and the insets in Fig.(b)-(d) show 
maps of the hole dispersion. The dark region is minimum of 
the dispersion. 



t" is sufficient to shift the LP to the disordered phase. 
The role of t' is less important. The shift of the LP is 
due to the tuning of the longer range hoppings. The 
"tuning" has been performed by nature in the cupratcs 
where the qualitative importance of t',t" is well known. 
These parameters give asymmetry between the hole and 
the electron doping. Holes go to the nodal points while 
electrons go to the antinodal ones resulting in dramat- 
ically different Fermi surfaces and magnetic properties. 
We follow nature and rely on the same mechanism. 

Our second major conclusion follows from the first 
one and it concerns SCS at the QCP. According to Fig 
OJa), at the QCP the lowest energy spectral function, 
k = (■§,§), does not have a pole, but only a cut. Ac- 
cording to previous studies for an immobile impurity, a 
cut indicates that the spin density is distributed in a 
power cloud around the hole [2(| [28| . Having a similar 
Green's function we directly project the results of Refs. 
2fil 28 1 to the present case. When approaching the QCP 



from the magnetically disordered phase the quasiparticle 
residue approaches zero, Z oc A 2 , z ps 0.4, as the triplon 
gap A approaches zero. The fraction of spin localized at 
the hole goes to zero oc Z. The rest of spin is distributed 
around the hole over a disk of radius i? oc 1/ A. At 7- -C i? 
the spin density is oc l/r Q , 1 < a < 1.5. Therefore the 
average radius of the spin cloud (r) ~ R diverges at the 
QCP, indicating SCS. On the other side of the QCP, deep 
inside the AF phase, the hole interaction with a magnetic 
field is described by pseudospin [H, 22|. This interaction 
implies a partial SCS [l7[. Evolution of the spin cloud 
when approaching the QCP from the AF phase is not 



clear at present. Because of the diverging magnon cloud 
the hole effective mass also diverges at the QCP. Draw- 
ing analogy with the cuprates we note that the effective 
mass measured in MQO (l5l | diverges on approaching the 
QCP identified by neutron scattering [18l420j . 

In conclusion. We consider the bilayer t — if — t" — J 
model with strong interlayer coupling J± . The hole dop- 
ing is low, so the magnetic dynamics are driven only by 
the interlayer coupling. At J±_ — > oo the hole Fermi sur- 
face is connected and centered at k = (tt, tt). At cer- 
tain J^_ p the AF correlations reconstruct the connected 
Fermi surface into four separate pockets centered close 
to k = (±^,±Tj). We have demonstrated that the LP 
where the topology of the FS is changed can be located 
within the magnetically disordered phase. The LP is 
driven by purely dynamic AF correlations in absence of 
any static magnetic order. The physics behind the LP is 
the following. Strong AF correlations make the nearest 
site hopping t almost idle, while the distant hoppings t' 
and t" are influenced to a much lesser extent. A balance 
of these effects gives rise to the LP. Because of the impor- 
tance of the AF correlations the LP is located not very 
far from the magnetic QCP. 

Wc have also demonstrated that if the LP is located in 
the magnetically disordered phase then, on approaching 
the magnetic QCP, the hole spin and charge separate. 
The separation scale is equal to the magnetic correlation 
length which diverges at the QCP. 

We are not aware of 2D materials with a low concentra- 
tion of charge carriers and with a magnetic QCP driven 
by a separate parameter. Perhaps such materials will be 
synthesised in the future and/or the model can be re- 
alised with cold atoms. However, the most important 
outcome of the analysis is the demonstration of the prin- 
cipal possibility of the topological transition and spin- 
charge separation. We believe that conceptually the re- 
sults are applicable to the single and multi-layer cuprates. 

We are grateful to A. V. Chubukov and G. Khaliullin 
for useful comments. Numerical calculations performed 
in this work were done at the National Computational 
Infrastructure at the ANU, under project u66. 
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